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Abstract The first stars are assumed to be predomi- 
nantly massive. Although, due to the low initial abun- 
dances of heavy elements the line-driven stellar winds 
are supposed to be inefficient in the first stars, these 
stars may loose a significant amount of their initial 
mass by other mechanisms. In this work, we study the 
evolution with a prescribed mass loss rate of very mas- 
sive, galactic and pregalactic. Population III stars, with 
initial metallicities Z = 10~^ and Z ~ 10^^, respec- 
tively, and initial masses 100, 120, 150, 200, and 250 Mq 
during the hydrogen and helium burning phases. The 
evolution of these stars depends on their initial mass, 
metallicity and the mass loss rate. Low metallicity stars 
are hotter, compact and luminous, and they are shifted 
to the blue upper part in the Hertzprung-Russell dia- 
gram. With mass loss these stars provide an efficient 
mixing of nucleosynthetic products, and depending on 
the He-core mass their final fate could be either pair- 
instability supernovae or energetic hypernovae. These 
stars contributed to the reionization of the universe and 
its enrichment with heavy elements, which influences 
the subsequent star formation properties. 

Keywords first stars, stars: models, evolution, mass 
loss 



1 Introduction 

It is widely accepted that the so called Population III 
stars, i.e. the first stars formed from the primordial 
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gas with an extremely low metallicity, have played an 
important role in the evolution of the universe. The 
first stars are assumed to be predominantly very mas- 
sive and they were responsible for reionization of the 
universe and its enrichment by heavy elements. 

The formation of the first stars and hence also their 
Initial Mass Function (IMF) is significantly influenced 
by the low metallicity of the protostellar material. Due 
to its low opacity, its gravitational attraction can not be 
effectively opposed by radiative pressure and the con- 
tinuing accretion gives rise to more massive stars. 

Numerical hydrodynamic simulations of star forma- 
tion generally confirm the tendency to higher masses of 
the first stars, although quantitatively the results differ 
co nsiderably. 



iBromm and Larson I ([200J) foun d the masses of the 



Omukai and Palla 



first s tars between 10^ and 10'^ Mq, 

( 2003[ ) followed the evolution of accreting protostars 
and found ^ 600 Mp) as the upper limit for massive 



stars. According to lAbel et al. l(j20Q2l) . the final masses 



are uncertain because a single molecular protostar seed 
of ^ ^Mq at the centre of a 100 core of a 
protogalaxy is formed. The cores may fragment into 
clusters of massiv e (MS) or very massive stars (VMS). 



Nakamura et al. I (jl999l ) found 3 Mq as the mass of 
the first stars which may grow up to ~ 16 Mq by the 
accretion. The IMF of Population III could thus be 
bimodal with peaks at ^ 1 — 2 Mq and ~ 10^ Mq 



Nakamura and Umemura I 120011 ). iTumlinson et al 



20041 ) defined strong VMS hypothesis ("The first gen- 



eration were exclusively VMS" ) and the weak VMS hy- 
pothesis ("The first generation included VMS in addi- 
tion to MS with M < 50 Mq"). 

The structure, evolution and properties of the first 
stars have been modelled i n several studies includin g 
the Paper I of this series ( Bahena and Klapp 2010l ). 



The results showed that the Population III stars are 
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smaller but more luminous and hotter than the normal 
metal-richer Population I an II stars of equal masses. 

The massive luminous normal stars have a high mass 
loss mainly due to the stell ar wind driven by ra diative 
pressure in lines of metals ( Castor et al. 19751 the so 
called CAK-wind). This process cannot take place in 



the metal- free stars ( Bromm and Loeb 1120031 ). This is 



why practically all models of first-stars evolution have 
neglected completely the mass loss. 

Several recent studies, however, revealed that even 
in the absence of CAK-wind there may occur a consid- 
erable mass loss due to different alternative processes. 

A reliable quantitative theory of the mass loss rate in 
the Population III stars is still missing, but it is of great 
interest to investigate the possible consequences of the 
mass loss on the properties of these stars. For this pur- 
pose we extend hereby our calculations of the conserva- 
tive case of the first-stars evolution in Paper I assuming 
a mass loss rate given by an ad-hoc parametrization. 

This work is organized as follows: In Sect. [2] the mass 
loss mechanisms are reviewed. In Sect.[3]we describe the 
initial conditions and results of the computed models, 
the main physical variables and their properties. Then, 
in Sect. m we discuss our results in comparison with the 
conservative case. Finally, in Sect. [5] we summarize the 
conclusions. 



2 Mass Loss Mechanisms 

The mass loss in stars of Population I and II has been 
studied both theoretically and observationally. How- 
ever, an extrapolation of the results to stars of Popu- 
lation III is questionable due to very different physical 
properties of both the interiors as well as the atmo- 
spheres of the stars. 

Regarding the short lifetimes of the massive first 
stars, their mass loss rates have thus to be estimated 
from the theoretical considerations and only an in- 
direct evidence can be gained from observations of 
later generations of stars. E.g., the abundances of 
eleme nts and isotopes observed in metal-poor halo 
stars (jBeers and Christlieb 1 120051 ) suggest that these 
stars were formed from an interstellar matter en- 
riched by mass lost during the H- and He- burn- 
ing in the first stars rather than by debris of pair- 
instability supernovae (SNe) exp ected to be the final 



stage of the first VMS evolut ion (jHirschi et al. 
Cescutti and Chiappini llioiol) . 



2008 



2.1 Stellar winds 

The mass loss via stellar winds influences significantly 
the evolution of massive Population I and II stars and 



contributes a great deal to the enrichment of interstellar 
matter by heavy elements. 

For cool winds of P opulation I red giants on the RGB 
Reimers I (jl975l . Il977l ) found an empirical relation giv- 
ing the mass loss rate as 



M 



L 

7r 



LR 



(1) 



where L, R, M and g are the lu minosity, radius, mass 



and s urface gravity of the star. ISchroder and Cuntz 
(|2005() obtained a semi-empirical generalization of the 
Reimers' relation. 

The stellar winds are significantly enhanced at the 
early-type stars by radiative pressure either in contin- 
uum or especially in lines. The dependence of the mass 
loss rate due to CAK-wind on the metallicity has been 
scaled by the relation 



(2) 



where the ex ponent varies between 0.5 — 0.8 for 0.01 < 
Z/Zp < 1 dLeitherer et al. Ill992l: [Vmk et al. Ill999 



200C 



2^ 1 . u i|ijeiuueiei eu ai. iiia a/iivmK. ni ai. i , 
a l200l[ ). However. iKudritzki I |2002l) has found that 



this relation breaks down below a certain threshold of 
th e metallicity. 



Krticka and Kubat 



( 


2006 




2009): 



Krticka et al. 



stars. They showed that weak stellar winds may be 
due to absorption in lines of CNO elements carried by 
mixing from the interiors to the atmospheres of the 
stars. 

2.2 Rotation 

The Population HI stars are supposed to be faster rota- 
tors than the stars with highe r metallicity ( Maeder et al. ^ 



19991: iMevnet et al~l l2006bl ). The num erical hydro 



dynam ic models of star-formation by IStacv et al. 



1 20111) also indicate high angular momenta of the first 
stars, what infiuences their structure, evolution and 
final fate. The effects of rotation on the structure 
and evolution o f massive stars were studied by the 
Geneva Group (IMevnet and Maeder I 119971: iMaeder 



1998; 
200C, 



Maeder and Mevnet 1 1200 It IMevnet and Maeder . 
20021 ). The rotation mav increase the mass loss 



by several processes: 

1. The rotational mixing enhances the atmospheres 
with heavier e lements which facilita te the radiatively 



driven winds ( Lamers et al. 1119951 ). 



The gravitational darkening enhances the radiative 
flux and the mass loss on poles of the rotationally 
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obla te stars while reducing the anKular-momentum 
loss (Maeder and Desiacques II2OOII : iDwarkadas and Ow( 



2002 ) 



The role of magnetic fields in activity of late-type 

C8ttkr^ and their winds is generally known. The grow- 



3. The centrifugal force reduces the gravity on the 
equator. When approaching the critical rotation, 
the decretion disks are formed through which the 
angula r momentum and significant amount of mass 
is lost (iKrticka et al. IboiOal boill) . 



2.3 Pulsation 

The pulsational instability of very massive stars might 
contribute somewhat to the total mass loss rate al- 
though it is not clear what the contribution will be. 
Cool supergiants do indeed have winds which are ini- 
tiated by pulsations, but then further support, such as 
radiative force on the dust grains, is needed to sustain 
a significant outflow. 

The stability of very massive stars wa s analyzed and 



compa red with metal enriched stars by iBaraffe et al. 
( 20011 ). Mass loss was not taken into account because 



they assumed it to be negligible for zero-metallicity 
stars. A linear non-adiabatic analysis was carried out. 
In their analysis, they considered that the stars stabi- 
lize as soon as they evolve from the ZAMS, then one 
can expect decreasing mass loss rates along the main 
sequence. This is, pulsations do not lead to significant 
mass loss on the main sequence. For their calculations 
they used the so called flux-freezing convective theory 
that assumes that perturbation in the luminosity is only 
given by the radiative luminosity. However, the stabil- 
ity properties of the models co uld change if convec tion 



20051) 



is properly taken into account (jKlapp et al. 

Thus, the issue of whether very massive stars could 
be vibrationally unstable to pulsation and induce pulsa- 
tionally driven mass loss is still open. The contribution 
of pulsation to mass loss remains uncertain. 

2.4 Other mechanisms 



Smith and Owocki I (j2006l ) suggested that the mass loss 
during the evolution of Population III stars may be 
dominated by optically thick, continuum-driven out- 
bursts or hydrodynamical explosions, similarly as in 
the very massive Luminous Blue Variable (LBV) stars 
and r] Car in particular. Unlike the steady line- 
dri yen winds, this mechanis m, investigated also, e.g., 
by van Marie et al. |2008al lbl). is insensitive to metal- 



licity but requires a high, nearly Eddington luminosity, 
close to which the mass loss rates of the wind highly in 



crease (Kuc 


ritzki 


Vink et al. 


I2OII) 



2002t iGrafener and Hamann 112002 



This condition is well satisfied by 
the first stars owing to their high masses and luminosi- 
ties. 



ing observational evidence of magnetic fields in early- 
type stars (e.g. Be-stars) suggests that they could also 
be present in the first stars and consequently infiuence 
thei r winds or rnass loss via their decretion disks (cf.. 



Puis et al. I boost lud-Doula et al. IboOSl and ref- 



e.g. 

erences therein). However, it is not clear whether there 
was a magnetic field in the interstellar medium in the 
early universe. 

Finally, because the first stars were also created in 
binaries and multiple systems, they could be subjected 
to mass-exchange, mass loss from their systems and en- 
hancement of their winds by tidal forces from compan- 
ion stars as it is theoretically well understood and ob- 
servationally confirmed for the present generations of 
stars. 

Motivated by these studies we explore a scenario for 
the evolution of very massive and low-metallicity Pop- 
ulation HI mass losing stars. 



3 Evolutionary Models at Low Metallicity 

3.1 Input physics and initial conditions 



In agreement with Castellani ( 2000l ) we take as the 
Population HI stars all the stellar objects with metal- 
licity below a not yet well defined upper limit between 
Z = 10~^° and 10~^. Depending on their metallicity 
we distinguish between galactic and pregalactic stars, 
the former with metallicity Z = 10~^, and the later 
with Z = 10~^, both considered as lower metallicity 
stars, but not zero-metallicity stars. 

The concept of critical metallicity has been used 
to characterize the transition between Popu lation HI 
and Population II star formation modes (lOmukai 



200lHBromm et al. Il200"li ISchneider et al 112002112003 



Mackev et al. 1 120031 : iBromm and Larson II2004I ) 



Our evolutionary models for both galactic and pre- 
galactic Population HI stars have been calculated 
by using a stellar code described in Paper I by 



Bahena and Klapp 



()2010l ). These models assume the 
standard quasistatic non-rotating configuration of the 
stars. A semiconvcction treatment for convective trans- 
port of energy has been used. In addition to the as- 
sumptions on the underlying physics used in Paper I, 
we include now also a mass loss parametrization. 

In view of the uncertainty in the dominant mass 
loss process in the Population HI stars, the corre- 
sponding mass loss rates and their changes with the 
stellar parameters as summarized in the previous Sec- 
tion, we have adopted a simple parametrization of the 
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m ass loss rate, based on the first relations proposed 



bv ILucv and Solomon I (jlQTOf ) and iBarlow and Cohen 
( 19771 ) to explain the mass loss for hot luminous stars, 



M = NL/c^ 



(3) 



where L is the luminosity of the star, c is the speed of 
light, and TV is a mass loss parameter. This parameter 
is taken as an average measure of the mass loss rate and 
it is kept constant in the course of the whole evolution. 
It is based on the simplifying assumption that the mass- 
energy output of the star in the form of non-zero rest- 
mass particles is proportional to its radiative output. 

Different models of stellar winds give the mass loss 
rates in the form of power-law in luminosity with var- 
ious exponents and they usually also depend on other 
stellar parameters. Let us note that unlike the Rcimers' 
formula ([1]) corresponding to the stellar winds of late- 
type stars, which for chosen M gives the M pro- 
portional to R increasing for orders of magnitude for 
stars evolving from the main sequence to the red-giant 
branch, our approximation ([3]) intended to an unspec- 
ified mass loss mechanism varies linearly with L only, 
i.e. much less in the course of evolution of the star. 

It can be expected that the mass loss rate is more 
complicated in reality, however, to get an insight into its 
possible effects on the stellar evolution before a proper 
model of the mass loss mechanism will be available we 
are limited to the use of such an ad hoc assumption. In 
our parametrization, the mass loss parameter N = 100 
gives for the luminosity L ~ IO^Lq typical for stars 
with mass of the order 100 Mq a mass loss rate of M 
lO~^M0yr^^ which is able to reduce significantly the 
initial mass of the star during its characteristic life- 
time r ^ lO^yr. Calculations with different N values 
should be indicative for the way in which the evolution 
of massive first stars is affected by the mass loss. 

To facilitate a comparison with models without mass 
loss the initial conditions for the models presented here 
have been chosen the same like in Paper I. 

The chemical composition corresponding to galactic 
lower metallicity stars is {X, Z = 0.765, 10^^), and for 
pregalactic ones {X, Z = 0.765, 10^^) by using a value 
Y = 0.235 as a primordial helium value. The initial 
stellar masses have been chosen 100, 120, 150, 200, and 
250 M0. For the mass loss parameter we have used the 
values TV = 50 and 100. 

3.2 Properties of the models 

The evolutionary models presented here predict the ef- 
fective temperature, luminosity and other physical vari- 
ables and quantities for very massive Population III 



mass losing stars during the hydrogen and helium burn- 
ing phases, for the initial masses and mass loss param- 
eters chosen. 

Figures [1] to |4] show some physical variables during 
the hydrogen and helium burning for pregalactic stars 
with masses 100, 120, 150, 200, and 250 Mq, metalhcity 
Z = 10~^ and mass loss parameter N = 100. 

Population III stars are hotter than their metal en- 
riched counterparts, and their locus in the HR-diagram 
is shifted to the left upper part; pregalactic stars are 
bluer than the galactic ones. 

For a given stellar mass of main-sequence stars its lo- 
cation in the HR-diagram depends strongly on metal- 
licity but also on the mass loss strength. During the 
hydrogen burning, the higher the mass loss rate, the 
greater the reduction in luminosity. 

Our models exhibit a high production of ionizing 
photons. This confirms the important cosmological 
consequences of the first generation of stars for the 
reionization of the universe. 

Depending on their initial mass and the mass loss 
rate the final fate of our models could be as supernovae 
or hypernovae explosions. 

a) Central density and temperature 

Lower metallicity stars are denser than the enriched 
stars. At the beginning of the evolution, galactic and 
pregalactic stars have different central densities which 
increase with decreasing Z. The central density in- 
creases with decreasing stellar mass and it is higher for 
lower metallicity stars. 

As the evolution goes on, the mean molecular weight 
increases and the density and temperature increase. 
During hydrogen burning, both the central density and 
temperature increase slowly until hydrogen is exhausted 
in the core. 

The most massive stars, and lower metallicity 
stars, have higher ce ntral temperature. According to 



Bromm et al. I (|2001l ) the central temperature is related 



to the mass and radius as 



R 



(4) 



For stars with M > 100 Mq the central temperature 
is insensitive to both mass and radius and depends upon 
the composition only. In the absence of metals, Tc is 
higher and the star self-generates the small amount of 
CNO ele ments needed for energy generation to support 
the star (ICastellani Il2000l) . 
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Fig. 1 Central density for 100 Mq (solid line), 120 Mq 
(dash-dot-dot-dot-dash), 150 Mq (dashes), 200 (dash- 
dot-dash) and 250 Mq (dots) pregalactic Population III 
stars with metallicity Z = 10^^ and mass loss parameter 
A'^ — 100 during the hydrogen burning phase. 
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Fig. 2 Idem dito. Luminosity 
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Fig. 3 Idem dito. Effective temperature 



b) Effective temperature and radius 

At the beginning of the main sequence, galactic stars 
with metallicity Z = 10~^, and mass loss parameter 
iV = 50, have an effective temperature T^s = 69582, 
73416, and 75586 K for a 100, 150 and 200 Mq, respec- 
tively. For the same stellar masses, pregalactic stars 
with Z = IQ-^ have Toft = 91865, 97109, and 100090 
K, respectively. Due to mass loss, the effective temper- 
atures decrease during core hydrogen burning. At the 
end, these temperatures are lower than for the conser- 
vative case. 

On the main sequence, 100, 150, and 200 Mq galac- 
tic stars with metallicity Z = 10^^, have initial radii 
R = 7.67, 9.54 and 11.15 i?©, respectively, and for pre- 
galactic stars with Z = 10~^, their radii arc R = 4.48, 
5.55 and 6.46 i?©, respectively. 

With decreasing effective temperature, the radius in- 
creases. For stars evolving with high mass loss param- 
eters, their radii at the end of the hydrogen burning 
increase even more. 

The changes of the radius and effective temperature 
can be understood in terms of changes in the convective 
core size. To satisfy the boundary condition 

L = TiacR^T^s (5) 

the effective temperature decreases with increasing ra- 
dius. An expanding convective core implies a decreas- 
ing radius and increasing effective temperature. 

h) Convective core size 

During the hydrogen burning the convective core 
evolves in different ways depending on whether the 
stars are mass losing or not. If a star evolves with mass 
loss its convective core decreases continuously during 
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Fig. 4 Idem dito. Radius 
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the hydrogen burning phase. In the conservative case, 
a hehum core is formed only at the end of hydrogen 
burning when the star is contracting during its transi- 
tion to hehum burning. 

At the end of hydrogen burning, the convective core 
size, defined as the ratio qcc ~ Mcc/M of the convective 
core mass Mcc to the total stellar mass M, is larger for 
evolution with mass loss than for the conservative case. 
This is because for mass losing stars their initial mass is 
reduced during the evolution. Figure[5] shows the evo- 
lution of the convective core size during the hydrogen 
burning for pregalactic Population III stars and mass 
loss parameter = 100. 

d) Luminosity and Gamma factors 

Luminosity increases or decreases during evolution with 
mass loss depending on the mass loss rate. The rate of 
decrease of the luminosity is higher for large mass loss 
parameters. 

In Figure [Bl the mass-luminosity (M — L) relation 
is shown for 100 Mq galactic and pregalactic Popula- 
tion III mass losing stars during the hydrogen burning. 
For low values of the mass loss parameter = 50, the 
luminosity increases during hydrogen burning but de- 
creases for large values of A^ = 100. This is because in 
VMS the opacity is mostly due to electron scattering. 
Then, as the evolution goes on, the opacity decreases 
and the mean molecular weight increases. Both of these 
effects contribute to an increase of luminosity which is 
proportional to the stellar mass. The M ~ L relation 
reads 
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Fig. 5 Convective core size for 100 MQ (solid line) , 
120 Mq (dash-dot-dot-dot-dash), 150 A/q (dashes), 200 Mq 
(dash-dot-dash) and 250 Mq (dots) pregalactic Population 
III stars with metallicity Z = 10~® and mass loss parameter 
N = 100 during the hydrogen burning. 
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where the value of the exponent 77 depends on the mass 
loss rate. The evolution reaches a given mass with 
higher luminosity, the lower the value of A^ is. 

Figure [7] shows the Gamma factor during the core 
hydrogen burning, for pregalactic Population III stars 
with a mass loss parameter A^ — 100. 

Galactic, and pregalactic Population III stars evolve 
during hydrogen burning near the Eddington upper lu- 
minosity limit. With and without mass loss, the ratio 
between the luminosity of the star and its Eddington 
luminosity is similar. 

In Figures |8] and 9, the Gamma factor is plotted 
as a function of the radius showing the phase when 
this factor approaches to unity, during the hydrogen 
burning, for mass losing galactic and pregalactic stars, 
respectively. 



Fig. 6 Mass-luminosity relation for 100 Mq galactic and 
pregalactic Population III stars with metallicity Z = 10~^ 
and Z — 10^^, respectively, and mass loss parameters N — 
50 and N — 100, during the hydrogen burning. 
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Fig. 7 The F factor for pregalactic Population III stars 
with metallicity Z = 10~^ and mass loss parameter N = 100 
during the hydrogen burning. 
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Fig. 8 The F factor as function of tiie radius for galactic 
Population III stars with metallicity Z = 10~® and mass 
loss parameter TV = 100 during the hydrogen burning. 
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Fig. 9 Idem dito. Pregalactic Population III stars with 
metallicity Z — 10~^ and mass loss parameter A'' = 100 
during the hydrogen burning. 
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Fig. 10 logT - logp plane for 100, 120, 150, 200 and 
250 Mq pregalactic Population III stars with metallicity 
Z = 10~^, and mass loss parameter A*" — 100, during the 
hydrogen and helium burning phases. 



e) p — T plane 

In the p^T plane the galactic and pregalactic very mas- 
sive Population III mass losing stars occupy the upper 
loci of a nondcgenerate and nonrclativistic gas. In this 
zone there is an upper boundary at which pair produc- 
tion effects could become important. The most massive 
stars tend to approach to it, however, our models are 
below this boundary. Regarding the gas in thermody- 
namic equilibrium, mass losing stars studied here are 
dominated by radiation pressure. Fig. [10] shows the 
evolution of conditions in the centres of stars. 

f) HR diagram 

As in the conservative case, the locus of very massive 
Population III stars in the HR diagram is in the left up- 
per part. Pregalactic stars are hotter than the galactic 
ones. The most massive stars are the most luminous. 
Depending upon the star's initial mass and mass loss 
rate, its evolution can be different. 

During hydrogen burning, the luminosity of a star 
without mass loss increases while a helium core is 
formed. Then, as the star evolves, the luminosity re- 
mains almost constant while the effective temperature 
decreases. However, for different stars their luminos- 
ity and effective temperature decrease with decreasing 
mass and increasing metallicity. The luminosity is lower 
for mass losing stars than for stars without mass loss. 

Different stars evolve with different lifetimes. De- 
pending on the mass loss rate, both the initial hydro- 
gen and helium burning phases take place at the blue 
side of the HR diagram. When the stars evolve with 
assumed mass loss rates they move toward the red, like 
in the conservative case. With and without mass loss 
the stars do not experience the AGB phase. 

Figures [Tl]to [T4]show evolutionary tracks in the HR 
diagram for 100, 120, 150, 200 and 250 M©, for mass 
loss parameters = 50 and A^ = 100, corresponding 
to galactic and pregalactic Population HI stars with 
metallicity Z = 10~^ and Z = 10^^, respectively. 



4 Discussion 

Population HI stars are defined as all metal deficient 
stars which share the common peculiarity of experi- 
encing self-producti on of C sometime s during their H- 



burning evolution ( Castellani 2000t) . That is, these 



stars may build-up a small fraction of C nuclei {Zq ~ 
10"^") vi a the triple-a reactions durin g the pre- main 



sequence ( Cassisi and Castellani 19931 ) 



Then, these stars produce their nuclear energy dur- 
ing the H-burning phase from a combination of the 
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Fig. 11 Evolutionary tracks in the HR-diagram for 100, 
120, 150, 200 and 250 Mq galactic Population III stars with 
metallicity Z = 10~^ and mass loss parameter A*' — 50, 
during the hydrogen and helium burning phases. 
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Fig. 12 Idem dito. Mass loss parameter A'' = 100. 
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Fig. 13 Evolutionary tracks in the HR-diagram for 100, 
120, 150, 200 and 250 Mq pregalactic Population III stars, 
metallicity Z = 10^^, and mass loss parameter A'' = 50, 
during the hydrogen and helium burning phases. 



6.8 


1 1 1 1 1 


6.6 


250 M 














^ 




200 M 

















^ 6.4 
cn 

O 


150 M 








6.2 


- 120 
100 M 






N=100 - 

z=io-^ 


6.0 


1 




1 


5 


1 4.8 


4.5 


4.2 3 



bg L„ 

Fig. 14 Idem dito. Mass loss parameter A'^ = 100. 
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proton-proton (pp) chain and the CNO-cycle with the 
small amount of C previously built. To provide the nec- 
essary luminosity, the star has to reach very high central 
temperatures ^ 10^ K for the simultaneous occur- 
rence of the helium burning via the Sa process. This 
burning produces certain amount of heavy elements to 
follow efficiently the CNO-cycle. 

In our models, the lower is the metallicity, the sooner 
is activated the triple-a reaction during the main- 
sequence. Lower-metallicity pregalactic stars without 
mass loss produce higher nuclear energy via the triplc-a 
reaction from the beginning of the main-sequence. For 
galactic stars contribution of this reaction is at the end 
of the main-sequence only. This is shown in Figure [T51 
For mass losing stars, the nuclear energy production 
via the pp-chain, CNO-cycle and 3a reaction is higher 
than in the conservative case, and 3a reaction starts 
from the beginning of the main-sequence, as it is shown 
in Figure [TBI 

As a result of lower energy-generation rate in the 
core, the stars maintain high central temperatures to 
support the stellar mass against the gravitational col- 
lapse. This high temperatures in the core and reduced 
opacity in their envelopes make these stars hotter and 
more compact than their enriched counterparts, de- 
pending on the initial stellar mass and metallicity. 

The lower is the metallicity, the more compact and 
hotter arc the stars. Table 1 shows models with 100 Mq 
when they settle down on the main-sequence. A galac- 
tic star with metallicity Z = 10^® has a radius of 
7.68 i?o, while a pregalactic star with Z = 10~^ has 
4.49 or with Z = 10"^° it has 3.69 i?©. Concern- 
ing the effective temperature this is 69582, 91865 and 
101800 K, respectively. On the ZAMS, the pregalactic 
lower metallicity stars are denser, hotter and slightly 
more luminous than galactic ones. 

With effective temperatures of Toff ^ 10^ K these 
stars are very efficient at producing p hotons capable 



of io nizing the hydrogen and helium (jBromm et al. 
200 ll) . 



Lower metallicity stars are shifted to the blue upper 
part of the HR-diagram. When stars are gradually en- 
riched with heavy elements they evolve redwards. This 
is shown in Figure [17] for a 100 Mq model without a 
mass loss. Mass losing stars evolve reducing their mass, 
luminosity and effective temperature, as it is shown in 
Figure [m 

In our models, evolutionary tracks describe different 
paths on the HR-diagram depending on stellar mass 
and mass loss rates. Galactic and pregalactic Popula- 
tion III mass losing stars move from the left to the right 
in the upper part of this diagram. With higher mass 
loss rates, the luminosity and the effective temperature 
decreases more than in the conservative case. 




Fig. 15 Nuclear energy generation via pp-chain, CNO- 
cycle and 3a-reaction for galactic and pregalactic Popula- 
tion III stars with initial metallicity Z = 10"*^ to Z = 10"^°, 
without mass loss during the hydrogen burning. 
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Fig. 16 Nuclear energy generation via pp-chain, CNO cy- 
cle and 3a-reaction for galactic and pregalactic Population 
III stars with initial metallicity Z = 10~^ and Z = 10~^ , 
with mass loss during the hydrogen burning. 



Table 1 Some characteristics of lower-metallicity 100 Mq 
stars when settle down on the main-sequence. Models be- 
come more hotter, denser, luminous and compact, as metal- 
licity is decreasing. 



z 


logTc 


log Pc 


, L 
log -j— 

Lq 


log Tcff 


log R 


10-" 


7.84650 


0.87274 


6.08715 


4.84250 


11.72763 


10-^ 


7.92150 


1.09791 


6.09132 


4.88074 


11.65323 


10-* 


8.00027 


1.33451 


6.09652 


4.92088 


11.57558 


10-8 


8.08323 


1.58381 


6.10314 


4.96315 


11.49433 


IQ-i" 


8.17087 


1.84730 


6.11174 


5.00775 


11.40942 
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As a result of their high mass and temperature, very 
massive stars are dominated by radiation pressure and 
they have luminosities closer to the Eddington limit. 
The specific luminosity LEdd ~ -W, and the total 

luminosity depends on the total mass of stars onl y but it 



is independent of the IMF (jBromm et al. 1120011) , which 



be a dominant effect in massive star evolution. They 
estimate that low Z stars may loose mass due to fast ro- 
tation as a consequence of their initial s low winds, and 



is a unique characteristic of very massive stars. 

Mass loss can affect the stellar evolution in a variety 
of ways. In the present work, nuclear lifetimes for mass 
losing stars with a mass loss parameter TV = 50 are 
similar to those in the conservative case. 

Hydrogen and helium burning lifetimes are decreas- 
ing with increasing stellar mass. On the other hand, 
the nuclear lifetimes are increasing with increasing mass 
loss rates. Nuclear lifetimes also are different depend- 
ing on the initial chemical composition of the stars. For 
M > 100 galactic and pregalactic Population III 
stars the nuclear lifetimes during the hydrogen burning 
are th < 3 x 10^ years. The helium burning lifetimes 
are shorter than the hydrogen burning lifetimes. 

As discussed in Section [51 it is still an open question 
if the lower metallicity stars could have lost a substan- 
tial part of their mass and which could be the mecha- 
nisms. We argue that a mass loss is necessary to have 
a good mixing and transport of the produced elements 
to the surface of the star. In the conservative case, 
without any mass loss, the mixing is very deficient and 
the CNO elements are not transported to the surface. 
However, in a scenario with mass loss, the mixing is 
very efficient. The main result is a production of N^"*, 
which originates primarily in the pregalactic star s, and 
subsequently in the galactic ones (jBahena 1 120061 ) . as it 
wi ll be shown in details in a fo rthc oming paper. 

Mevnet and Maeder ( 2002 ) and Maeder and Mevnet ^ 
( 20031 ) lave shown that the stellar rotation and mass 



initia l rotation velocities may be faster (jMaeder et al 
1999t) . Their models of massive stars at low Z do not 



reach in general the red stage. However, rotation per- 
mits the envelope to largely expand and the star moves 
to the red supergiant stage, modifying all the mod- 
els outputs. Then, the rotation affects the evolution 



of very metal-poor massive stars (iMevnet et al. II2005I) 



and ca n induce mass loss in many ways (j Mevnet et al. 



2006a|). The compactness of low metallicity stars is a 



loss can affect the predictions of the stellar yields espe- 
cially for He, C, N and O. At low Z, the rotation may 



key property which has many effects in relation with 
stellar rotation. 

For the present work. Tables 2 and 3 indicate some 
features for galactic and pregalactic Population HI mass 
losing stars. The first column lists the initial stellar 
mass Mi, the second one refers to the mass loss pa- 
rameter N. The next columns indicate the lifetime th 
during hydrogen burning given in Mega- years, the mass 
ejecta during the hydrogen burning AMh, the convec- 
tive core size 5cc,h during this phase and the corre- 
sponding He-core mass Mne. Then there follow the 
lifetime tho during helium burning in Mega-years, the 
convective core size gcc.Hc during this burning phase 
and the CO-core mass Mco- The mass ejecta, the He- 
and CO-core masses are given in solar masses. 

In our models the mass ejecta are important during 
hydrogen burning and they increase with increasing ini- 
tial stellar mass. For 100 and 200 Mq pregalactic stars, 
with N = 100, the ejecta AMh are 26.82 and 59.08 Mq, 
respectively. 

In the present work, mass losing stars form a convec- 
tive core size slightly greater than in the conservative 
case. Stars evolving with mass loss parameter N = 100 
form a helium core Mro = 32.98, and 71.27 Mq for 100 
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Fig. 17 Evolutionary tracks in the HR-diagram for a 
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Table 2 Nuclear life times and masses during hydrogen and helium burning for 100, 120, 150, 200 and 250 Mq galactic 
Population III stars with initial metallicity Z = 10"*^, and diverse mass loss parameters. 
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26 82291 


0.45071 


73 17709 


32 98165 


3 1 5890 


30761 




22 48512 


120 





2.89075 


0.00000 


0.43025 


120.00000 


51.64200 


2.93901 


0.40415 


120.00000 


48.49800 




50 


2.87198 


17.78249 


0.47813 


102.21751 


48.87326 


2.87981 


0.33142 


102.21992 


33.87773 




100 


2.94026 


32.84446 


0.47004 


87.15554 


40.96659 


2.94880 


0.34399 


87.04112 


29.94127 


150 





2.67959 


0.00000 


0.44097 


150.00000 


66.14550 


2.72591 


0.42386 


150.00000 


63.57900 




50 


2.67557 


23.56230 


0.49297 


126.43770 


62.32999 


2.68253 


0.38855 


126.42807 


49.12363 




100 


2.73644 


43.49372 


0.49057 


106.50628 


52.50440 


2.74280 


0.38651 


106.33505 


41.09956 


200 





2.46549 


0.00000 


0.45237 


200.00000 


90.47400 


2.50979 


0.42152 


200.00000 


84.30400 




50 


2.43380 


32.71329 


0.50490 


167.28671 


84.46306 


2.48056 


0.37267 


167.27864 


62.33973 




100 


2.51988 


59.07631 


0.50571 


140.92369 


71.26652 


2.52582 


0.34120 


140.68414 


48.00143 


250 





2.32990 


0.00000 


0.46402 


250.00000 


116.00500 


2.37348 


0.45782 


250.00000 


114.45500 




50 


2.34590 


41.44418 


0.51272 


208.55582 


106.00500 


2.35321 


0.40578 


208.37607 


84.55484 




100 


2.38727 


45.41756 


0.51194 


174.58244 


89.37573 


2.39382 


0.41550 


174.27963 


72.41319 



and 200 M© galactic {Z = 10^^) stars, respectively. 
For pregalactic (Z = 10^^) stars, A/hc = 32.29 and 
68.72 Mq, respectively. The corresponding carbon core 
masses are 22.48 and 48.00 M© for the first case, and 
32.13 and 70.24 M© for the second one, respectively. 

Studies of core-collapse SNe have discovered two 
distinct types of supernovae explosions: (i) very en- 
ergetic SNe or hypernov ae (HNe), and ( ii) ve ry faint 
and low energy SNe. In iNomoto et al."l ( 2004 ) it was 
proposed that the first generation supernovae were 
the explosions of ~ 20 — 30 M© stars and some of 
them produced C-rich, Fc poor ejecta. The ejecta 
of these explosions can well account for the abun- 
dance pattern of EMP stars. In contrast, the ob- 
served abundance patterns cannot be explained by 
the explosions of more massive, 130 — 300 M© stars. 
These stars undergo pair-in stability SNe explosioris and 



are disrupted completely (jUmeda and Nomoto 
Heger and Wooslev l[20oi ). 



2002 



If all stars evolve with mass loss they can reduce 
their initial mass to massive stars range, as our models 
show. Depending of their initial mass and mass loss 
rate, its final fate could be either SNe or HNe. The 
HNe explosions connected to some low-^ gamma-ray 
bursts coul d be relevant to early nucleo synthesis at low 
metallicity (jUmeda and Nomoto II2002I) . 



5 Conclusions 

The main aim of the present study was to construct 
models for very massive low-metallicity both galactic 



and pregalactic Population HI mass losing stars. We 
have calculated and discussed here the main properties 
of these stars in comparison with the conservative case 
presented in Paper I. Our results show that the mass 
loss could be very important for massive low metallic- 
ity stars, although the physics and mechanisms of the 
mass loss are not well known yet. According to the 
line-driven wind theory, the mass loss rate is negligible 
for the very low metallicity. However, the mass loss 
may occur due to other processes which arc metallic- 
ity independent, such as effects of continuum radiation, 
pulsation and/or rotation etc. In our approach we have 
chosen an ad hoc rate of the mass loss and we used 
the know physics to explore its effects on the stellar 
evolution at very low Z . 

For evolution with mass loss, the reduction of the 
stellar mass makes the central temperature to increase 
less rapidly than for constant-mass evolution. Thus the 
mass of the convective core decreases more rapidly as 
the evolution proceeds. However, the core mass frac- 
tion is larger in a star evolving with mass loss. The 
luminosity of a star evolving with mass loss is lower 
than for constant-mass evolution and in consequence 
its main-sequence lifetime is somewhat increased. 

The evolution of VMS both with and without mass 
loss takes place in the upper part of the HR-diagram. In 
the absence of heavy CNO elements, lower metallicity 
galactic and pregalactic stars contract and settle down 
on the ZAMS with central temperatures high enough 
for some '^He to be converted to CNO nuclides through 
the triple-a reaction and CNO cycles. The energy gen- 
erated stops the contraction and the star re-expands. 
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Table 3 Nuclear life times and masses during hydrogen and helium burning for 100, 120, 150, 200 and 250 Mq pregalactic 
Population III stars with initial metallicity Z = 10"", and diverse mass loss parameters. 
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0.00000 
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2.73923 


0.40604 


120.00000 
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50 


2.70386 


17.72365 


0.46846 


102.27635 


47.91238 


2.74578 


0.46763 


102.26871 


47.82392 




100 


2.76936 


32.99806 


0.46448 


87.00194 


40.41066 


2.81274 


0.46295 


86.43422 


40.01472 


150 





2.50358 


0.00000 


0.39613 


150.00000 


59.41950 


2.54207 


0.41884 


150.00000 


62.82600 




50 


2.52162 


22.89308 


0.48177 


127.10692 


61.23630 


2.52587 


0.49071 


127.09890 


62.36826 




100 


2.57387 


41.94500 


0.47805 


108.05500 


51.65569 


2.57775 


0.49219 


108.04709 


53.17970 


200 





2.30257 


0.00000 


0.41782 


200.00000 


83.56400 


2.34043 


0.43612 


200.00000 


87.22400 




50 


2.34086 


31.96341 


0.48552 


168.03659 


81.58513 


2.34750 


0.48872 


168.02637 


82.11785 




100 


2.38205 


58.43802 


0.48546 


141.56198 


68.72268 


2.38542 


0.49668 


141.41945 


70.24021 


250 





2.17765 


0.00000 


0.42354 


250.00000 


105.88500 


2.21375 


0.44310 


250.00000 


110.77500 




50 


2.22612 


40.76518 


0.48886 


209.23482 


102.28653 


2.23255 


0.49352 


209.16652 


103.22786 




100 


2.25313 


74.29857 


0.48770 


175.70143 


85.68959 


2.25955 


0.49271 


175.39363 


86.41820 



In all cases, lower metallicity pushes the stars to 
a smaller radius at higher density and temperature. 
Therefore, galactic and pregalactic Population III stars 
are more compact, hotter and denser than their metal 
enriched counterparts. Due to their different chemical 
composition, the luminosity is slightly lower for pre- 
galactic stars. 

The evolution of both galactic and pregalactic Popu- 
lation III stars is influenced by the increase of the mean 
molecular weight in the convective core and the effect 
of mass loss. The mean molecular weight affects the 
star by increasing their luminosity. As a result of the 
increased mean molecular weight in the convective core 
the central temperature also increases. With increased 
molecular weight the convective core shrinks. The con- 
traction of the convective core makes the star to expand 
and the effective temperature to decrease. The propor- 
tionality between mass and luminosity causes that the 
lu minosity decreases as the st ar lo ses mass. 



Heger and Wooslev I (|2002[ ) and lHeger et al. I (|2002 



200 



^ have studied the final fate of VMS considering 



only the conservative evolution. Stars with 140 Mq < 
M < 260 M© explode as Pair- Instability Supernovae 
Explosions (PISNe) causing complete disruption. In 
these cases, Mhc ~ 64 — 133 Mq. 

Our models for 200 Mq galactic and pregalactic stars 
develop in the conservative case Mne ~ 90 and 84 M©, 
respectively. With moderate mass loss rate (iV = 50), 
they develop Mne ^ 84, and 82 M©. With increasing 
mass loss rate (A'' = 100), stars develop Mhc = 71 and 
69 Mq. Then, these galactic and pregalactic stars could 
explode as PISNe. 



However, a c cording to scena r ios by Umeda and Nompto h 



2OO2I I2OO3I) : lUmeda et al~l |2002l ): iNomoto et al 



20041) ■ the stars with M < 130 Af© could explode like 
HNe. Our 100 Mq galactic and pregalactic models de- 
velop A^Hc ^ 42 and 38 Mq in the conservative case 
{N = 0), A/hc ^ 39 and 38 A/© with a moderate mass 
loss rate {N = 50), and with N = 100, they have 
Az/hc = 33 and 32 Af©, respectively. Then, models for 
< 100 Mq mass losing stars, which evolve with high 
mass loss rates, can reduce appreciably their initial 
mass during the hydrogen, and helium burning. Their 
cores could also explode like HNe, which arc likely to 
make an important contribution to the early galactic 
chemical evolution. 

The most important effect in the transition in the 
star formation properties from the first stars to present 
day is the presence of a trace amount of heavy ele- 
ments. The concept of critical metallicity has been 
used to characterize the transition between P opulation 
III and Population II star formation mo des ('Omukai 



2OOII: iBromm et al. Il200lt iBromm and Larson ..200 



f 



This important parameter value of critical threshold for 
enabling a formation of stars with lower mass is esti- 
mated to be Z ~ 10-3-5^© - 6.32 x 10"^. In our 
very massive evolutionary models, their initial mass is 
reduced when they lost mass reaching the range of mas- 
sive stars. Due to the mass loss an efficient mixing of 
CNO elements is reached. After hydrogen burning, the 
initial metallicity increases providing an initial enrich- 
ment of the IGM with heavy elements, giving place to 
the next generation of stars with increasing metallic- 



ity. In successive generations the EMP stars observed 
in present days could be formed. 

The evolution of the first stars with mass loss is still 
an open problem. It would be interesting to suppose 
that the rotation or other mechanisms driving the mass 
loss is simultaneously responsible for a mixing of ele- 
ments in the stellar interiors. The main problem, how- 
ever, is to find a realistic model for the mass loss rate 
due to the different mechanisms and their combinations 
and to include it into self-consistent models of the stel- 
lar evolution. 

This work has been partially supported by the Mexican 
Consejo Nacional de Ciencia y Tecnologi'a (CONACyT). DB 
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